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ABSTRACT 

In this work we show the results obtained from the VLBA data at 15 GHz 
of OJ287 in the super-resolution mode. The data showed a jet configuration 
in the form of a "fork" where superluminal components emerge via stationary 
components at the northwest and the southeast close to the core to form parallel 
trajectories along the southwest direction in the plane of sky. This agrees with a 
source structure of an extended, broad morphology of OJ287. 

Subject headings: Galaxies: Jets — BL Lacertae objects: individual (OJ287) — 
radio continuum: galaxies — techniques: interferometric 



1. Introduction 



OJ287 is a well-studied BL Lac object and distinguished in the optic al by its series of 
prorainent fiares (double p eaked) which recur with a period of ~ 12 years (jSillanpaa et al. 



19881 : IVillforth et al. II2010I ). OJ287 is the only known extragalactic source showing evidence 
of a major periodic component in the optical and bringing into a scenario of a binary black 
hole model. In radio, VLBI observations of OJ287 reveal a core-jet structure with a size not 
much greater than 1.5 mas. In 1995-96 the jet-position-angle is ~ -90° which rotates in a few 
years to southwest and showing a jet position angle of ~ -120° in 1998. The change of the jet 
direction was modeled as a ballistic jet prece ssi on taking into account the o ptical periodicity 
of 12 years ( Tateyama and Kingham |2004[) . Va.ltonen and Wiik I (12012 ) proposed a new 



model rescaling the nodding model of iKatz I (Il997f ) with a precession period of 120 years. In 
2008, O J287 presented a dramatic change on the structure of the source. VLBA observations 
at 43 GHz of lAgudo et al. I ( 120101 ) showed a strong, bright knot in the southeast of the core 
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in 2008, almost perpendicular to the previous dire ction of the jet. An analysis including 



earlier VLBA data at 43 GHz, lAgudo et al. I (|2012[ ) found that such transverse components 



were already present on the source which was interpreted as the crossing of the jet from 
one side of the line of sight to the other of the inner jet. In this work we investigate the 
struct ure of the inner j et of OJ287 at 15 GHz from 1995 to 2012 in the super-resolution 
mode (ITateyama 1120091 ). 



2. Data and Imaging procedure 



We have used the visibility VLBA data at 15 GHz from MOJAVE (Monitoring of jets 
in active galactic nuclei with VLBA experiments) program and 43 GHz from BU (Boston 
University). Table 1 list the VLBA data used in this work. The columns of the table are: 
1 - epoch, 2 - frequency, 3 - peak brightness, 4 - synthesized beam size, and 5 - position 



angle of the beam. The images were obtained using the Difmap Package ( IShephard 1119971 ) 
and processed in an automate mode using the mapping script of Taylor and Shephard. The 
restoring beam for both data was 0.1 mas. The pixel size of the imag es was 0.05 mas. Th e 
major features of the maps at 43 GHz have already been presented by lAgudo et al. I ( 120121 ). 



The CLEAN is a method of deconvolution invented by iHogbonj (119741 ). The main as- 
pects o f the CLEAN can be foun d in a confront with ME (Maximun Entropy) deconvolution 



m e.g., iRastorgueva et al.l ( 1201 ll ). One of the aspects of the CLEAN is that it shows better 
images with compact features as in the case of OJ287. The CLEAN algorithm produces 
a set of point components which convolved with the CLEAN beam and superposed on the 
residual map generates the CLEAN map. The synthesized beam (dirty beam) is obtained by 
averaging the outputs of all pairs of an array. The CLEAN beam is usually a fit of the main 
lobe of the dirty beam a nearly Gaussian with angular resolution k, \/ D (diffraction- limited 
telescope). Deconvolution methods dealing with structure details smaller than the standard 
beam (diffraction pattern) are so -called the super-resolution techniques. The interferomet- 
ric data have been always subject of deconvolution techniques to obt ain sharper stru ctures. 
In a weak way we can mention the super- resolution effects of ME M (IBaikova Il200ll ) or the 
inadequacy of the standard beam to image interferometric data (IReid II2006I ). In a strong 
way (large fact o r of s harpening) we can mention the direct study of the fringe visibility by 
Kovalev et al. I (120051 ) which have shown compact components r nuch smaller than standard 
beam and the super-resolution study at 15 GHz (factor 5) by ITateyama I (120091 ) showing 
structures comparable to the structures given by the resolution at 43 GHz. Such interfer- 
ometric data along with a set of maps containing adequate time sampling observations in 
order to recognize features of the source in adjacent maps produce a robust data results. 



- 3 - 



The CLEAN-point image was also examined by carrying out simulation of VLBA image 
using the MPS tasks UVCON and IMMOD. The UVCON generates visibilty data from a 
VLBA array geometry and the IMMOD modifies the existing map image by the addition of 
models. The simulated data were imaged using the Difmap Package. 



3. The structure of OJ287 at 15 GHz 

Figure 1 shows the VLBA observations of OJ287 in the super-resolution mode from 1995 
to 2012 at 15 GHz. The solid line on the maps is a reference line at position angle of -120° 
to guide the structure of the jet in the maps. We follow the identification of components 
adopted by Agudo et al. (2012) at 43 GHz. The bright jet feature at the eastern end of 
super-resolution maps denominated "C" in our maps corr esponds to th e core but it is not 



always coincident with position of the core proposed by lAgudo et al. I (120121 ) at 43 GHz 
maps. Accordingly the identification of component "a" which is very close to the core along 
southeast-northwest direction also differs. In addition to the components "C" and "a" a 
discrepancy on the component "d" was removed by introducing a feature "dl" in our maps. 
Also, using the 43 GHz data, new features "q" and "r" were included on the maps of 2011 
and 2012. 

In the 90's the source shows a narrow jet morphology dominated by two main stationary 
components. One stationary component is at ~ 0.3 mas denominated "A" and the other 
at ~ 1 mas denominated "G". The compone nt "G" can be id e ntifie d with superluminal 



component C6 on the VLBI Observations by iTateyama et al. I (1l999l ) before becoming a 
stationary component. The others superluminal components "K","0", "R","T", ... becomes 
prominent knots when passing at the location of the stationary feature at 1 mas. 

It is well known that a super-resolution deconvolution is very sensitive to produce ar- 
tificial structures. While a super-resolution image of a set of compact components in a 
rectilinear direction is easy to capture, the super-resolution image of a complex structure 
shows complicated artifacts as arcs. An interesting aspect of such artifacts is that they con- 
tain information on the source structure and this can be recovered performing a simulation 
of VLBA observation. The simulation can iteratively fit the structure of the source (a model) 
with observed images (super- resolution images). Figure 2 shows the model of the source for 
epochs 2005, 2009 and 2010. The feature on the map of 2005 Jul showing an arc blending 
with a stationary component at 0.3 mas is generated by collinear components along the jet 
(southeast) with a northern knot close to the core. The transversal component in 2009 Aug 
is formed with a knot in the southeast of the core and the transversal structure in 2010 May 
is produced by a southeast /core/northwest source configuration. 
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A remarkable feature on the super-resolution maps is the appearance of a strong, trans- 
verse component along the core in 2008-9. In fact it corresponds to the striking change on 
the structure of the source with appearance of a bright knot in the southeast of the core 
reported by Agudo et al. (2010, 2012). In addition to this bright transversal component in 
2008 other weaker transversal features can be recognized on the maps. The arc (the artificial 
feature of the super-resolution mode) on the core of 2002 Jun is related to the knot "L" on 
the 43 GHz maps of Agudo et al. (2012), in 2003a May is related with knot "S", in 2004 
Apr-Aug-Sep is related to knot "a" and in 2005 Jul is related to knot "a" in the north. 
These transversal features appear on the northwest or north of the core. In 2006 Apr and 
2008 May component "a " appears on the south of core and in 2010 May the knots "n" and 
"1" appear on the north of the core. Several knots especially early components as "L" , "S" 
or even "a" appear as bright feature but does not show superluminal motion, suggesting a 
stage of carving a new channel on the jet. 

In 2007 the transversal component seen at 43 GHz is not seen at 15 GHz. The component 
is very weak at 15 GHz or the resolution at 15 GHz is not sufficient to produce the feature 
in the southeast-northwest direction. But in May 2008 the transversal component appears 
on the southeast of the core and gradually becomes very prominent component at 15 GHz 
and 43 GHz. 



4. Position of the core 

The configuration of the structure at 15 GHz, in particular the transverse feature in 
2008-2009, indicates the component "C" as the proper position of the core. If we choose 
the component in the southeast of the core as the core it would deform the structure of 
the jet. We would have to presume that the position of the core has changed (using the 
optically thin jet components as reference). The structure of the CLEAN-point maps on the 
2004-2005 also indicates that the bright knots near the core appear in the north/northwest 
an d in the south /south east of the core and not from the southeast to northwest as adopted 



bv lAgudo et al. I (]2012f ). 



In this way we propose that the stationary components (bends on the fiow) are formed 
in the southeast and northwest near the core. The first knots would show bright feature near 
the core without producing superluminal components. The absence of superluminal motion 
would reflect the stage of carving a stationary component or distorting the flux of the jet 
to form a broad jet. On this interpretation the bright knots in the southeast (south) and 
northwest (north) of the core would show a form of a "fork" with a southeast /core/northwest 
conflguration where superluminal components would emerge. 
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Figure 3 shows the 43 GHz maps at resolution of 0.1 mas. The sohd hne on the maps is 
the same reference hne drawn on Figure 1 to guide the structure of the source and the origin 
of the curve is placed on the core accordingly with SE/core/NE configuration. The inner 
jet already appears in the southeast-northwest direction. In 2007, the core is the brightest 
component and the transversal component is weak or very close to the core. In 2008, the 
bright transversal knot is in the south of the core. 

The core on the configuration of the structure of the jet around 2000 at 15 GHz is the 
brightest component at the base of the line extending from the stationary components at 
0.3 mas and 1.0 mas. This configuration starts to change in 2001 with the formation of new 
stationary components in the southeast-northwest direction. The source evolves to show 
parallel trajectories with superluminal components emerging via stationary components in 
the southeast and northwest close to the core (bends on the flux flow). This agrees with a 
source structure of an extended, broad morphology of OJ287, produced by parallel trajec- 
tories than a broad transversal expansion of a single jet. Figure 4 shows the superposition 
of maps from 2008 to 2012 at 15 GHz and 43 GHz to emphasize the "fork" aspect of the 
source. 

The southeast /core/northwest configuration remov es the abrupt variation of amplitude 



of core in mid-2004 (see Figure 3 of lAgudo et al. II2012I ). On the scenario of a broad jet the 



trajectory of components is consistent with a prominent jet in the southwest direction. The 
jet-position-angle would show large dispersion only on the inner jet. Furthermore the core 
(in the sub-parsec broad line region) would be the site of the variable gamma-ray emission. 

The striking similarity of VLBA structure at 15 GHz and 43 GHz in 2011-2012 supports 
the southeast /core/northwest configuration rather placing the core on the southernmost 
component of the transversal core region. This implies that the location of the mm core is 
not different of the cm core. 



5. Conclusion 

In this work we show VLBA data at 15 GHz (MOJAVE) of OJ287 by examining the 
structure in the maps in a super-resolution mode along with 43 GHz VLBA maps. The result 
showed that the CLEAN-point maps at 15 GHz follow very closely the structure of the 43 
GHz maps. The innermost jet shows a configuration of SE/core/NW in the form of "fork". 
New components emerge from the core via components at southeast and northwest near the 
core to show superluminal components in parallel trajectories in the southwest direction in 
the plane of sky. 
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Fig. 1. — 15 GHz VLBA images of OJ287. The peak flux densities of the maps are indicated 
on the left side of the maps. Contour levels are 1, 2, 4, 8, 16, 32, 64, 128, 256 and 512 x 
0.003 Jy beam-^ 
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Fig. 1. — (Continued) 
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Fig. 2. — Simulation and Model maps at 15 GHz. Contour levels are 1, 2, 4, 8, 16, 32, 64, 
128, 256 and 512 x 0.003 Jy beam"!. 
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Fig. 3. — 43 GHz VLBA images of OJ287. The two strongest flux densities along the 
southeast-northwest direction are shown on the left side of the maps, the flux on the lower 
position corresponds to a knot in the southeast of the core. Contour levels are 1, 2, 4, 8, 16, 
32, 64, 128, 256 and 512 x 0.005 Jy beam-^ 
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V 




Fig. 4. — Superposition of maps at 15 GHz and 43 GHz to emphasize the "fork" aspect of 
the source. Contour levels at 15 GHz are 1, 2, 4, 8, 16, 32, 64, 128, 256 and 512 x 0.0045 
Jy beam-i and contour levels at 43 GHz are 1, 2, 4, 8, 16, 32, 64, 128, 256 and 512 x 0.005 
Jy beam~^ 
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Table 1. VLB A data at 15 GHz (MOJAVE) and 43 GHz (BU) used in this work. 
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Table 1 — Continued 
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Table 1 — Continued 
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-9 


2010 May 24 


15 


4.1 


1.3 X 0.6 


-12 


2010 Jul 09 


15 


4.3 


1.1 X 0.5 


-11 


2010 Aug 28 


15 


2.5 


1.0 X 0.5 


-5 


2010 Oct 18 


15 


4.5 


1.0x0.5 


-5 


2011 Feb 27 


15 


5.5 


1.0x0.5 


-3 


2011 Aug 26 


15 


5.0 


1.1 X 0.6 


-7 


2011 Dec 29 


15 


4.8 


1.1 X 0.6 


2 


2012 may 24 


15 


4.7 


1.0x0.5 
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Table 1 — Continued 



Epoch 


Frequency 


^peak 


Beam 


p. A. 




GHz 


T 1 —1 

Jy beam 


mas X mas 




Ofin? Tun 14 


4"^ 


1 


5 -X- 2 


-1 5 


2007 Till 12 


4"^ 


Q 


4 -X- 2 


-8 

(J 


2007 Aiip- flfi 


4"^ 


8 


5 -X- 2 


-17 


2nn7 Aiip- 


4S 


1 2 


4 x 2 

W . ^ -A. W . Zj 


-9 


9007 Spd 9Q 


4'^ 


1 




-1 9 




4'^ 


2 


5x09 


-91 

Z( i- 


9008 Tan 1 7 


4'^ 


2 


4x09 

U . t: -A. W . Zi 


-9 


2008 Feb 2Q 


43 


1.4 


0.4 X 0.2 


-16 


2008 Tun 12 


43 


2 7 


5x02 


-14 


2008 Till Ofi 


43 


2 2 


4 -X- 2 


-1 1 


2008 Aiip- Ifi 


43 


2 3 

Zj . O 


5x02 

L/. O -A. \j .£j 


-22 


2008 Sen 1 


43 


2 1 


5x02 

\J .XJ -A. \J .£j 


-1 6 


2008 Nov Ifi 


43 


2 Q 


4 -x- 2 

W . t: -A. \j .^j 


-9 


2008 Dpp 21 


43 


2 7 


4x02 


_7 

1 


200Q Tan 24 


43 




4x02 

W . t: -A. W . Z( 


-1 3 


200Q Fpb 22 


43 


4 7 


4x02 


-1 5 


200Q Anr 01 






U.t: -A U .Zi 


14 


200Q Mav "^0 


43 


3 5 


4x02 


_7 


200Q Tun 21 

£j\j\jXj U Llil Zi± 


43 


2 7 


3x02 


_i 


200Q Till 2fi 


43 


3 7 


3x02 


_i 


zUU9 Aug lo 


A O 


4.9 


U.4 X U.z 


■1 

-1 


2009 Sep 16 


43 


5.1 


0.4 x 0.2 


1 


2009 Oct 14 


43 


5.6 


0.4 x 0.2 


-6 


2009 Oct 16 


43 


7.4 


0.5 X 0.2 


-21 


2009 Oct 20 


43 


7.4 


0.5 X 0.2 


-16 


2009 Oct 25 


43 


5.4 


0.4 X 0.2 


-5 


2009 Nov 28 


43 


7.1 


0.5 X 0.2 


28 


2010 Jan 10 


43 


8.2 


0.4 X 0.2 


-4 


2010 Feb 11 


43 


3.3 


0.4 X 0.2 


-7 


2010 Mar 06 


43 


6.0 


0.4 X 0.2 


-4 
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Table 1 — Continued 



rijpocii 


rrec[uency 


T 

'-peak 


ijeam 


P A 




vjjnz 


jy Deani 


inao X niab 





2010 May 19 


43 


1.8 


0.4 X 0.2 


-11 


2010 Jun 14 


43 


1.8 


0.4 X 0.2 


-7 


2010 Aug 01 


43 


1.4 


0.3 X 0.2 


-7 


2010 Aug 21 


43 


1.4 


0.4 X 0.2 


-10 


2010 Sep 18 


43 


3.6 


0.4 X 0.2 


-16 


2010 Oct 24 


43 


.0 


0.5 X 0.2 


-15 


2010 Nov 01 


43 


3.4 


0.5 X 0.2 


-4 


2010 Nov 06 


43 


3.6 


0.5 X 0.2 


-8 


2010 Nov 13 


43 


2.8 


0.4 X 0.2 


-5 


2010 Dec 04 


43 


3.1 


0.4 X 0.2 


-6 


2011 Jan 02 


43 


3.4 


0.5 X 0.2 


-17 


2011 Feb 04 


43 


3.6 


0.4 X 0.2 


-14 


2011 Mar 01 


43 


3.7 


0.4 X 0.2 


-12 


2011 Apr 21 


43 


3.7 


0.3x0.2 


3 


2011 May 22 


43 


4.9 


0.4 X 0.2 


-16 


2011 Jun 12 


43 


2.9 


0.3 X 0.2 


-2 


2011 Jul 21 


43 


2.7 


0.4 X 0.2 


-6 


2011 Aug 23 


43 


2.9 


0.5 X 0.3 


-18 


2011 Sep 16 


43 


2.1 


0.4 X 0.2 


-9 


2011 Oct 16 


43 


3.6 


0.4 X 0.2 


-5 


2011 Dec 02 


43 


2.7 


0.4 X 0.2 


-8 


2012 Jan 27 


43 


3.2 


0.3 X 0.2 


-7 


2012 Mar 05 


43 


3.4 


0.4 X 0.2 


-10 


2012 Apr 02 


43 


3.0 


0.3 X 0.2 


1 



